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We propose and analyze a mesoscopic device producing on-demand entangled pairs of electrons. 
The system consists of two capacitively coupled Mach-Zehnder interferometers implemented in a 
quantum Hall structure. A pair of electron wave-packets is injected into the chiral edge states of 
two (of the four) incoming arms; scattering on the incoming interferometers splits the wave-packets 
into four components of which two interact. The resulting interaction phase associated with this 
component leads to the entanglement of the state; the latter is scattered at the outgoing beam splitter 
and analyzed in a Bell violation test measuring the presence of particles in the four outgoing leads. 
We study the two-particle case and determine the conditions to reach and observe full entanglement. 
We extend our two-particle analysis to include the underlying Fermi seas in the quantum Hall device; 
the change in shape of the wave-function, the generation of electron-hole pairs in the interaction 
regime, and a time delay between the pulses all reduce the degree of visible entanglement and the 
violation of the Bell inequality, effects which we analyze quantitatively. We determine the device 
settings optimizing the entanglement and the Bell test and find that violation is still possible in the 
presence of the Fermi seas, with a maximal Bell parameter reaching B = 2.18 > 2 in our setup. 

PACS numbers: 73.23.-b, 03.67.Bg, 85.35. Ds, 73.43.Lp 



I. INTRODUCTION 

Quantum entanglement is a genuine property of quan- 
tum mechanics that has attracted a lot of attention re- 
cently due to its potential usefulness as a computational 
resource. As was first noticed by Bell' , entangled states 
can violate a certain type of inequality expressed in terms 
of the correlation functions of measured outcomes. Later, 
Clauser and coworkers- suggested a more transparent in- 
equality which was experimentally violated by photonic 
entangled states . 

The controlled creation and manipulation of electronic 
entangled states in solid state devices is a challenging 
problem and no experimental demonstration of a Bell in- 
equality violation with electrons is available so far. The 
main challenge lies in the fact that, in contrast to pho- 
tons, electrons are charged particles and thus strongly 
interact with the electromagnetic environment, leading 
to a fast decay of the coherence of the entangled state. 
On the other hand, the Coulomb interaction allows one 
to easily create entanglement among electrons. 

During the last decade, several strategies have been 
suggested to create entangled states of ballistically prop- 
agating electrons in mesoscopic devices. The initial pro- 
posal was to use the generic spin-singlet entanglement 
of Cooper pairs in a superconductor''', where the con- 
stituent electrons of a Cooper pair are injected into two 
different normal leads where they can propagate and 
thereby separate ballistically. Recently, the adiabatic 
splitting of Cooper pairs into normal conductors has 
been demonstrated experimentally . The constituents of 
Cooper-pairs are entangled both in the orbital (energies) 
and spin degrees of freedom. The advantage of the spin- 
entanglement over the orbital entanglement is that the 



spin coherence time can approach about 1 ms in semi- 
conductor devices. However, the reliable demonstration 
of spin entanglement between electrons in a Bell test' 
requires the detection of spin polarized currents at ar- 
bitrary polarization angles, that is still beyond the level 
of present technology. Alternatively, one can use the en- 
ergy degrees of freedom of the injected Cooper pairs as 
an entangled variable" . The corresponding Bell test then 
requires the measurement of time-resolved current corre- 
lators on the GHz scale in the frequency domain. 

Other proposals for two-electron entanglement make 
use of interacting quantum dot systems in the Coulomb 
blockade regime" or ballistic electrons in integer quantum 
Hall devices'"'' ' . The latter proposals involve the entan- 
glement of the electrons' orbital degrees of freedom and 
the corresponding Bell test requires the measurement of 
spinless current correlators among different leads of the 
setup. However, in contrast to the setups in Refs. 4,-^),!), 
in these schemes the entanglement is produced due to 
a postselection at the moment of measurement and no 
interaction among electrons is required, see Refs. 12 for 
the discussion of the role of projection. In the postse- 
lection schemes, entangled states are produced with less 
then 100 % efficiency, in contrast to the proposals using 
interactions, where the entangled states are prepared de- 
terministically by the unitary evolution. 

Another disadvantage in the above proposals is that 
they operate at finite bias voltage where Cooper pairs or 
electrons are injected stochastically into the device. As a 
result the Bell test requires the measurement of current 
correlators at short times (of the order of the voltage 
time oc h/eV). To go beyond the short time limit, one 
alternatively can operate in the tunneling regime where 
electrons are injected rarely and one can measure correla- 
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tors at larger times oc h/ {TeV) (T is the transparency of 
the injecting lead). However, actual applications rely on 
the controlled creation of entangled states on demand. It 
was suggested in Ref. ' '! to use Lorentzian voltage pulses 
to inject electron pairs on demand; the Bell test then 
requires only the detection of the total number of the 
transmitted particles through the different leads of the 
setup, and no short-time analysis is needed. The con- 
trolled injection of individual electrons into a device has 
been achieved recently, see Ref. 14. 

Considerable progress has been made in the fabri- 
cation of electronic Mach-Zehner interferometers where 
electronic transport occurs through the edge states of an 
integer quantum Hall state. Experiments^ have un- 
covered an Aharonov-Bohm interference pattern at high 
visibility ^ 62 %. Later, the interference of two indepen- 
dent electrons coming from different sources has been re- 
ported' ' , reproducing the original Hanbury Brown and 
Twiss experiment in optics with a visibility approaching 
70%. 

Several theoretical proposals have been made to use 
coupled electronic Mach-Zehnder interferometers to en- 
tangle two^'^''^'' or three independent electrons"'". All 
these proposal deal with an idealized situation where only 
two or three electrons are present in the setup. However 
in real experiments^ a non-monotonic dependence of 
the visibility factor has been observed at high bias volt- 
ages, that cannot be explained within a framework of 
non-interacting electron transport. Thus, the electron 
interaction cannot be ignored and indeed its proper ac- 
counting- ^ is required in order to explain the experimen- 
tal data. 

The purpose of the present work is to analyse an en- 
tangling device for (two) electrons consisting of two cou- 
pled Mach-Zehnder interferometers, where electrons are 
injected on demand on top of the filled Fermi sea. In 
contrast to earlier proposals here the injected elec- 
trons are interacting not only with each other but also 
with the underlying Fermi sea, leading to a parasitic en- 
tanglement with the surrounding electronic environment. 
We formulate the corresponding Bell inequality test and 
find the dependence of the Bell parameter ;B as a func- 
tion of the Coulomb interaction strength, the energy of 
the injected electrons, and the ratio between the mu- 
tual and self-capacitances of the interacting leads. As 
expected, the interaction with the Fermi sea leads to a 
strong decoherence and renders the Bell inequality vio- 
lation more challenging (but still doable) as compared 
with the idealized two-electron situation: injecting nar- 
row (high-energy) wave packets and for a mutual capac- 
itance surpassing the self-capacitance of the interacting 
leads we find that a Bell parameter B ~ 2.18 > 2 violat- 
ing the Bell inequality B <2 can be reached. 

In the following, we first discuss our setup and de- 
fine our goals, the calculation of the maximally possible 
Bell parameter quantifying the entanglement in the wave 
function and the implementation of the Bell test with the 
expected outcomes for B. In Section HI we discuss the 




FIG. 1: Two Mach-Zehnder interferometers capacitively cou- 
pled via the arms 1 and 2. Lorentzian-shaped pulses are in- 
jected into the left arms, entangled through the Coulomb in- 
teraction, and analyzed in a Bell test on the outgoing chan- 
nels. In spite of the decoherence generated by the presence 
of the Fermi sea, the Bell inequality can be violated, demon- 
strating the possibility to generate 'useful' entanglement cor- 
relations in this setup. 

Bell test with only two particles in the setup. The anal- 
ysis is extended to include the underlying Fermi seas in 
the device in Section IV and we conclude in Sec. V. Tech- 
nically, our analysis makes use of the fact that the effects 
of capacitive interaction between chiral electrons can be 
described by the action of voltage pulses-- . 

II. BELL SETUP 

We consider two Mach-Zehnder interferometers, the 
upper with leads '1' and '3' and a lower one with leads 
'2' and '4', where the two adjacent arms '1' and '2' are 
capacitively coupled, see Fig. 1. Each interferometer con- 
sists of two non-refiecting incoming and outgoing beam 
splitters A {A') and B {B') for the upper (lower) inter- 
ferometer. The beam splitters A and B {A' and B') can 
be described by the unitary transfer matrices t{a),t{/3) 
(t(a'),t(/3')) in the lead basis {'!', '3'} ({'2', '4'}), with 

x/ N COS a —sin a 
t{a) = . (1) 
^ ' sm a cos a 

parameterized by the angle a and corresponding expres- 
sions for the other matrices. 

We describe the interaction between the electrons 
propagating in the adjacent leads '1' and '2' by the 
Hamiltonian 

Hint ^ I J2 E,,N,N,, (2) 

ije{l,2} 

where Afi — J dx Ki{x) ■.'^\{x)i! i{x) : is the excess elec- 
tron number in a finite region of lead i defined by the 
coordinate kernel Ki{x) localized near the origin. We 
also use the notation : A:— A — {vac\A\vac) with \vac) 
denoting the ground state of the system. The coupling 
matrix E describes the coupling constants due to the self- 
capacitances of the interferometer arms (diagonal terms 
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Ell and £^22) and due to the mutual capacitance between 
the adjacent arms of the different interferometers (off- 
diagonal terms E12 = £21). For simplicity, we assume no 
interaction in the other leads '3' and '4', E33 = E44 = 0. 

The excess electron charges are injected into the in- 
terferometers in a well defined non-entangled quantum 
state. Due to the particle interaction, the outgoing state 
becomes orbitally entangled (in the particle number mea- 
sured in the outgoing leads) and this entanglement can be 
detected by performing a Bell test and observing the de- 
gree of violation of the corresponding CHSH-inequality. 
In the following, we formulate the Bell inequality in 
terms of cross-correlators (NiNj) {i G {'1','3'} and 
j e {'2', '4'}) between the number of excess particles 
transmitted through the different outgoing leads of the 
setup, Ni = J dt ■.'^l{x,t)'^i{x,t):. These correlators can 
be tuned by changing the magnetic fluxes $ and pen- 
etrating each interferometer as well as the angles /3 and 
/?' of the outgoing beam splitters (we will use symmet- 
ric splitters with a = a' = 7r/4 in the incoming leads). 
Defining the Bell correlator 

i.,,,(a>,c>0^«^^-^-)(^--^-)^ (3) 

{{Ni + Ns){N2 + N4)) 
and the corresponding Bell observable, 

S = i?^^,($,<i>')-fi?0^, ($,$') + 

+E^f,,{^,^')-E^^,i^,^'), (4) 

one can formulate the two-party CHSH-inequality, 



161 < 2. 



(5) 



The violation of the above inequality for some specific 
values of the magnetic fluxes $ = {$,$',$,$'} and an- 
gles /? = {f3, 13' , P'} quantifies the entanglement shared 
between the two interferometers. We will show that vio- 
lation (although not maximal in the full setup including 
the Fermi surfaces) can be attained by changing only the 
magnetic fluxes at fixed (optimal) values of the beam 
splitter angles, (3 — f3 and = (3' , reducing the num- 
ber of adjustable parameters, that might be helpful in a 
realistic experiment. 

The actual value of the Bell parameter depends on the 
settings of the Bell test, the magnetic fluxes <& and the 
adjustment of the outgoing beam splitters deflned by the 
angles (3: B ^ B{^,f3). In order to flnd the optimal 
setting of the Bell test where the Bell parameter attains 
its maxima, one needs to maximize B{^,/3) as a func- 
tion of all tuning parameters in the test. It has been 
shown in Ref. ' that this optimization problem can be 
done analytically by calculating the two-particle density 
matrix of the entangled particles. We deflne the two- 
particle density matrix for our setup by placing the out- 
going beam splitters B and B' into the asymptotic re- 
gion where ki.2{x) — > with no interaction present (as 
described by the limit x — > -l-oo). In this region (before 



the second beam splitters B and B'), the (many-particle) 
entangled scattering state \BB') freely propagates in the 
positive direction and we deflne the (reduced) density 
matrix 



[pBB'ln'n' OC 



dxdy{BB'\^{x)^,iy) (6) 
xi>,,iy)^,{x)\BB'), 



where (x) denote the free single-particle electronic held 
operators in the internal arms of the interferometers with 
i,j e {'1', '3'} and G {'2', '4'}. According to Ref. 23 
the maximal value of the Bell parameter is given by, 



BmaxipBE') = 2-\/ Ai + A2, 



(7) 



where Ai and A2 are the two largest eigenvalues of the 
3x3 symmetric matrix Tj^^^Tp^^, \ with [fp]nm = Tr {p- 
((Tn ® cfm)} and where (t„, n,m G {x, y, z} are the Pauli 
matrices. 

The phase accumulation due to magnetic fluxes $' 
and the scattering on the outgoing beam splitters act as 
a unitary rotation of the 4x4 density matrix. 



Pbb' 



Pout = ?7($, /?, /?') PBB' UH'^>, /?, /?') (8) 



in the lead (or pseudo-spin) basis. This unitary rota- 
tion does not change the entanglement of the state since 
it involves only independent single-particle rotations in 
each interferometer, ?/($, /3, /3') = ?/($, /3)(^J7($', /3') 
where. 



C/(<i>,/3) = £(/?) 



gi*/2 









(9) 



Below, we will use Eq. (7) in order to flnd the maximal 
possible value of the Bell parameter in the device and 
Eq. (8) in order to calculate the Bell parameter (4) from 
the correlators (3). 



III. BELL TEST WITH TWO ELECTRONS 

We consider first the idealized situation where only 
two single-electron wave-packets are injected into the in- 
coming leads '1' and '2' and no other electrons present 
in the system. The operator = J dx f{x)^l^{x) cre- 
ates a single-particle state with wave-function f{x) in 
lead a. Then, the incoming state with two wave-packets 
f{x) and g{x) in the leads '1' and '2' has the form, 
|in) = flg2 \vac). After scattering on the incoming beam 
splitters A and A', the wave-packets are split between 
the different internal arms of the interferometer. Before 
the wave-packets reach the interaction region the state of 
the system is a non-entangled product state. 



\AA') = (sin a/] + cos a fl] 



X [sm a + cos a §2 



(10) 
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Depending on the path chosen by the electrons, the 
various components of this product state evolve differ- 
ently and the overall two-particle state becomes non- 
separable with respect to the different interferometers. 
For the empty vacuum state the interaction effects ap- 
pear only in the scattering component where the elec- 
trons propagate through the adjacent arms '1' and '2'. 
Before the second beam splitters B and B' , the two- 
particle state has the form, 

\BB') = [cos a cos q;'(S' /I g|) -)- sin a sin a'(/|^4) 
+ sin a cos a' {flgl) + cos a sin a' {flgl)] \vac), (11) 

where S is the evolution operator corresponding to the in- 
teraction Hamiltonian (2). In order to find the scattering 
wave-function Sflg^lvac) of the two electrons propagat- 
ing in the interacting leads '1' and '2', we determine the 
eigenvalues and eigenvectors (ei,ei) and (£2,62) of the 
coupling matrix E in the interaction Hamiltonian (2) 
(note that the self-interaction terms have to be set to 
zero in the two-particle problem discussed here. En — 
for I = 1,2). Introducing the electron number operators 
■hi = X]j=i 2 ^ij-^j i^ij the j-th component of the vec- 
tor Bi) one can rewrite the interaction Hamiltonian in the 
quadratic form, 



£2 



(12) 



Next, we perform a Hubbard-Stratonovich transforma- 
tion to decouple the charge operators fii (although not 
really necessary here, this procedure is convenient when 
dealing with the many electron case later on, see Ap- 
pendix). The evolution operator then can be written in 
the form. 



S — DziDz2 exp 



dt{eizl{t) + e2zi{t)) (13) 



i / dt{eizi{t)ni{t) + e2Z2{t)h2{t)) 



xTexp 



where zi,2(t) are two real auxiliary fields obeying Gaus- 
sian statistics with {zi{t)zj{t')) = {i/ei)6ijS{t — t'). This 
transformed evolution operator describes the indepen- 
dent propagation of two electrons in different leads, 
each subjected to a time-dependent scattering potential 
V^{x,t) = V,{t) K,(x), i = 1,2, with 



Vi(t) = eizi(t) eu + f-2Z2{t) e2i- 



(14) 



Having mapped the two-particle evolution to two inde- 
pendent single-particle problems, we then solve the cor- 
responding Schrodinger equations for the chiral electron 
modes with linear spectrum e = vhk. The resulting two- 
particle wave-function '^i2{x, y;t) = {x,y \ S flgl\vac) 
emerging behind the interaction region takes the form 



i{x,y;t) = f{xt)g{yt)(ex:p\-i j dt' (15) 
X [Vi{t')Ki{xt + vt') + V2{t')n2{yt + i'i')] }), 



where Xt = x—vt and j/t = y—vt axe ballistically retarded 
variables, and the average has to be taken with respect 
to the Gaussian fields Zi{t). Carrying out the averaging 
of the exponentials, we find the result 



^I2{x,y;t)^f{x,)g{yt)e-^''^^(^^y\ 



(16) 



^12ix,y) = E12 J dT Ki{Xr)K,2{yT)- (17) 


The phase $12(0;, y) describes both the energy exchange 
and the deformation of the wave-functions f{x) and g{y) 
due to the interaction between the adjacent arms. In 
the end, the asymptotic {t — > 00) form of the interacting 
component in the scattered state (11) is given by the 
expression 



Sflgl^ I dxdy f{x)g{y)e 



-i*i2(2;,a).f,t / 



^\{x)^l(y). (18) 



Next, we calculate the two-particle density matrix (6) 
for the scattering state (11) inside the interferometer, 
PBB' =T:^x,v{\BB'){BB'\], 



PBB' = (Sq ® Sa') 



- 1 


1 


1 


V- 


1 


1 


1 


V 


1 


1 


1 


V 


_ V* 


V* 


V* 


1 . 



{Sc 



(19) 



with ~ diag{sinQ;,cosa} and 

V = I dxdy|/(a;)n5(y)|2e-*-(-'«). (20) 



We choose a particular form of the interaction kernels, 
Ki{x) = K2{x) = exp(— |a;|/a), where 2a is the length of 
the interaction region. Then the asymptotic form of the 
phase $12(2;, y) at x, y — -l-oo is given by. 



<l>,2(x,,) = ^oexp(-^)(l + ^), (21) 



with 



-El 2 To 



(22) 



the interaction-induced phase accumulated by the par- 
ticles during the simultaneous propagation through the 
interaction region; here, tq = a/v denotes the ballistic 
travelling time through the interaction region. We con- 
sider the case of simultaneous injection and choose in- 
coming wave-packets of Lorentzian form with ^ > 0, 



fix) = g{x) = \/ I 



TT X + 



(23) 



Then the parameter V can be expressed as a function 
of two dimensionless parameters, the phase <poi see Eq. 
(22), measuring the strength of the Coulomb interaction. 
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and 7 = C/fi quantifying the width of the incoming wave- 
packet, 



V(<^o,7) 



dx exp[^iipoe-^'>\='\ (1 + 27|a;| 

TT x'^ + 1 



(24) 



Consider first the situation with infinitely narrow in- 
coming wave-packets (7 — > 0) simuhaneously injected 
into the device. The relevant values of the phase $12(2;, y) 
in V are those near x ~ y, $12(2;, x) = ipQ, resulting in 



V(i^o,0) exp(-ji^o), 



(25) 



see Eqs. (20) and (21) (when the wave-packets are in- 
jected with a time delay Td, the phase ipc ^ ^uivr^) ^ 
ipo exp(— wTd/a) is reduced). The density matrix (19) cor- 
responds to a pure state with p^^, — pbb'- According 
to Eq. (7), the maximal value of the Bell parameter is 
given by, 



2 W 1 + sin^ {2a) sin^ {2a') sin' 



2 V'o 



(26) 



The maximal violation Smax = 2^/2 is attained for sym- 
metric incoming beam splitters with a = a' = 7r/4 and 
</?o = TT, requiring proper tuning of the interaction do- 
main. 

Next, we determine the optimal configurations of mag- 
netic fluxes and angles (<lJ,$,/?,^) and ($', /?', ^') 
for the upper and lower interferometer where the Bell 
inequality is maximally violated (we consider only the 
setup with symmetric incoming beam splitters). To do 
so, we derive the explicit form of the Bell correlation func- 
tion E'^^/ (<i>, $'), see Eq. (3). The correlation functions 
{NiNi') entering into E are given by the diagonal ele- 
ments of the two-particle density matrix pout after scat- 
tering at the outgoing beam splitters, see Eq. (8). Then 
the Bell correlation function (3) assumes the form (we 
can set (3 = /3 and /?' = (3' and still violate the Bell 
inequality maximally) 

^ sin 213 cos 2/3' [V cos{^ + (pc) - cos $] 

^ cos 2/3 sin 2/3' [V cos($' -I- .^c) - cos $'] 

I- sin 213 sin 2/3' [V cos($ + $' + (^c) +cos($- 



E 



(27) 



-*')], 



where V — \V\ and ipc — — argV. For infinitely narrow 
wave-packets simultaneously injected into the device, we 
have V — 1 and (pc — ifio, see Eq. (25). The Bell param- 
eter then takes the form 

B ^2 [sin 2/3 cos 2/3' sin $ + cos 2/3 sin 2/3' sin $'] sin y 

-I- sin 2/3 sin 2/3' [cos $ cos $' + cos $ cos $' (28) 
-I- cos $ cos !>' — cos !> cos $'] , 

where we have shifted all fluxes by (pc/2, $ — > $ + (^c/2 
and so on. The maximal violation Eq. (26) is reached 
when 
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FIG. 2: Bell parameters Smax and Bo as a function of of width 
parameter 7 = ^/0 at (/Jo = 7r for an optimal setting of Bell 
parameters (Smax, solid line) as given by Eq. (7) and for the 
non-optimal settings in Eq. (29) {Bo, dashed line). 



and with angles /3 = tt/4 and /3' ~ arcctg[sin((^c/2)]/2 for 
the outgoing bean splitters. Note, that one can approach 
the maximal violation fixing the angles /3 and /3' of the 
beam splitters to the above values and changing only the 
magnetic fluxes. 

In a next step, we consider the situation where the 
incoming wave-packets have a finite width and thus a 
finite 7. Here, in contrast to the ^ — case, the in- 
teraction deforms the initial shape of the single parti- 
cle wave-functions. This deformation appears only when 
both electrons are transmitted through the interacting 
arms '1' and '2' and is absent in the other scattering chan- 
nels. Thus the orbital entanglement in the scattering (or 
pseudo-spin) degrees of freedom is recorded in the specific 
form of the orbital wave-functions. Since our detection 
scheme captures only the total number of the transmit- 
ted particles and is insensitive to the specific shape of the 
wave-packets, the reduced density matrix pbb' in Eq. 
(19) is mixed and the pseudo-spin entanglement is re- 
duced (for the case of narrow wave-packets with ^ — 0, 
the tracing over x and y generates a mixed state as well, 
but does not reduce the entanglement which is only in 
the lead-indices) . A further reduction in entanglement of 
extended wave packets is due to the electrons not meeting 
each other in the finite-range interaction region, thereby 
reducing their effective interaction. 

When discussing the quantitative result for the Bell 
test with finite-width wave packets, we first analyze the 
value of the Bell parameter for the settings (29) which 
have been optimized for the 7 = case, 



Bo = {l + V)Jl 



2 '^c 



(30) 



($,!>) ^ (V2,0), ($',$') ^ (0,7r) 



(29) 



In the most favorable situation where ipc = tt, the Bell 
inequality cannot be violated for V < V2 — 1. In reality, 
the interaction phase ipc monotonically decreases from tt 
as 7 increases; in the limiting case of spatially extended 
wave-packets 7 — >■ 00, the electrons are unlikely to meet 
each other in the interaction region and ipc ^ 0. 

The settings in Eq. (29), however, are not optimal 
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FIG. 3: Contour plot of the Bell parameter Smax for optimal 
settings of the Bell parameters as a function of the interaction 
parameter ipo (horizontal axis) and the width parameter 7 
(vertical axis). For small- width, high-energy wave packets 
(7 0) the Bell inequality is maximally violated for odd 
multiples of n. 



when 7 is finite. In Fig. 2, we compare the maximal 
possible value of the Bell parameter B^ax deriving from 
Eq. (7) with the non-optimal value Bq as given by Eq. 
(30) choosing an interaction parameter ipo = tt. We find 
that Smax > Bq and the Bell inequality can be violated 
as long as 7 < 7c « 0.87 (the critical value for Bo is 
7c ~ 0.497). 

The maximal value of the Bell parameter in the entire 
parameter space ((/3o,7) is shown in Fig. 3. The region 
where the Bell inequality is violated {B > 2) is shrinking 
as the strength of the Coulomb interaction (as encoded 
in ipo) grows. Indeed, for a finite 7, a stronger inter- 
action implies a larger deformation of the wave-packets 
and thus more pseudo-spin entanglement is recorded in 
the shape of the wave-packets. As a result, the detectable 
pseudo-spin entanglement is reduced, thus lowering the 
Bell inequality violation. On the contrary, for 7 = 0, the 
maximal Bell parameter is a periodic function of ipo with 
maxima at ipo — tt, 37r, . . . , see Eq. (26). 



IV. BELL EXPERIMENT WITH FERMI SEA 

We now extend our analysis to the more realistic situ- 
ation including the (zero temperature) Fermi sea. First, 
we study the scattering problem where a single-particle 
Lorentzian wave-packet crosses an interaction region in a 
single lead. This problem has been studied before in Ref. 
22 and we restate the results here. Next, we consider 
the scattering of two Lorentzian wave-packets propagat- 
ing in two adjacent interacting leads and make use of 
these results to find a many-electron scattering state in 
the complete system with two coupled Mach-Zehnder in- 
terferometers. Knowledge of this state then permits us to 
determine the level of Bell inequality violation reachable 



in the presence of a Fermi sea. 



A. Scattering of a single wave-packet 

Consider a single-electron state with a wave-function 
f{x) created at — ^ —00 in the interaction-free asymp- 
totic region of a one-dimensional conductor. 



I/) = / dx fix)i'\x,to)\^t 



(31) 



where |$f) is the zero-temperature Fermi sea. We as- 
sume that all Fourier components of f{x) — J^k fk^^''^ 
are vanishing below the Fermi momentum, fk — for 
k < kp. For the specific case of a Lorentzian wave packet, 
see Eq. (23), the incoming state |/) can be created by ap- 
plying a local voltage pulse of Lorentzian shape. 



V{t) 



2h 



1 



er l + (i-to)7T2' 



(32) 



in the asymptotic region xq — > —00 of the conductor; 
here t — ^/vp is the duration of the pulse'^''. Then the 

state I/) can be written as^"', 



\f) ^um]F^\^T. 



(33) 



where is an electron ladder operator increasing the 
number of electrons in the system by one. The unitary 
operator tJ[(j>] describes the evolution of the wave func- 
tion under the action of the Lorentzian voltage pulse (32) 
applied at the position x — xq, 



U[(j)] = T+ exp ivp / dt (f>{t)p{xo,t) 



(34) 



with T-|_ the forward time-ordering operator, p{x, t) ~ 
: 'ii^ {x,t)^{x,t) : the electron density operator, and 
(j){t) — J dt'eV{t')/h is the phase accumulated by the 
electrons. 

The wave packet then propagates ballistically towards 
the interaction region around a; = and arrives there 
at t = 0, assuming xq = vptQ. Due to the interaction, 
the excess electron can exchange energy with the Fermi 
sea, thereby exciting electron-hole pairs. Subsequently, 
this many-particle scattering state (the excess electron 
plus the electron-hole cloud) ballistically propagates to 
the interaction-free region at large x — > +00 where it 
assumes the form 



I/) =r+exp 



dtHintit) I/) 



(35) 



with Hint given in Eq. (2). 

At first glance the scattering state |/) should be a 
rather complicated entangled state involving an excess 
electron and additional electron-hole pairs. However, the 
central result of Ref. -'2 says that the state |/) can again 
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be produced by applying an additional voltage pulse 
to the incoming state |/), |/) = C/[x]|/)i with 
X{t) ^ y dt' eV^{t')/h. The remarkable consequence of 
this fact is that the scattering state has the form of a sim- 
ple Slater determinant and is in fact non-entangled, since 
it can be produced by an evolution operator generated by 
a single particle Hamiltonian. 

The peculiarity of the Lorentzian incoming state be- 
comes clear when one goes to the bosonic formulation, 
expressing the fermionic field operator through bosonic 
creation/annihilation operators h\ and 6^, ^(x) oc 
exp[i ^t^yi^ibke-'^'^^ +b\e~^''''') / \/k]F , see Appendix for de- 
tails. Integrating over the coordinate in the definition of 
the incoming state (see Eq. (31)) with the Lorentzian 
form of ,f{x), the wave-packet can be written in the form 



|/)cxexp(-*^^fel)Ft|<i> 



fc>0 



(36) 



where |$b) is the bosonic vacuum state. This state 
can be recognized as a coherent state of bosons |/) 



nfe>o bfe) with amplitudes Vk 



l^fk. In the 



bosonized picture any electron interaction of the form 
V{x,y)p{x)p{y) is quadratic in bosonic variables and 
thus corresponds to a potential scattering of the bosons. 
Due to the chiral nature of the present scattering prob- 
lem, the back reflection of bosons is forbidden and the 
only result of that scattering is the appearance of a 
momentum dependent forward scattering phase dsc{k), 
bk — >■ fefc exp[j(5sc(fc)], where the actual form of dsc{k) de- 
pends on the particular form of the interaction kernel 
V{x,y). Hence, the scattered state of bosons is again a 
coherent state, 



I/) oc exp|^-z^ 



fc>0 



,-feC-ll5ec(fe) 

^/k 



(37) 



Going back to the electronic picture, this scattered state 
can again be created by applying a corresponding voltage 
pulse V{t) given by 



B. Scattering of two wave-packets 

Let us next analyse the scattering problem for two elec- 
tron wave-packets propagating in the two adjacent inter- 
acting leads ('1' and '2') of the upper and lower interfer- 
ometer (unprimed and primed) . We assume that initially 
at to ^ the incoming state is given by. 



1/5) = / dxdx'f{x)g{x')^l{x,toM{x',to)\^F). (39) 



In order to generalize the result of the previous section 
to the case of the two-particle scattering problem, we 
calculate the overlap of the scattered state \fg) = S\fg) 
and the state \fg)u (for later convenience we use x ^^nd 
x' referring to the upper and lower interferometers rather 
than the lead indices '1' and '2'), 



\f9)u^Ui[x]U2[x']\f9), 



(40) 



obtained by applying two distinct voltage pulses V^{t) 
and V^'{t) in the corresponding (upper and lower) leads 
to the incoming state \fg) in the non-interacting system. 
The overlap, 

OfMt),x'{t)] = {fg\ul[xMx']S\fg), (41) 

is a function of the phases x{t) ^-nd x'{t) ^-nd we want to 
find its maxima, Ofg = max,,,,^' [©[x, x']| as a function 
of phases x{t) and x'{t)- 

As shown in the Appendix, there exist particular 
phases x(0 ^-nd x'{t) producing an overlap (41) of unit 
magnitude, provided that the incoming wave-packets 
f{x) and g{x) have a Lorentzian form (the more gen- 
eral statement actually is that one can reach full overlap 
for any voltage-pulse generated incoming states). Thus 
we conclude that the scattered state \fg) — S\fg) can be 
obtained by applying two voltage pulses to the incoming 
state of the non-interacting problem. 



V(t) ~ - / duj cxp\~TUj + i6sciuj/vF) + iiot]+C.c. (38) 
e J 



Even more, for an initial voltage pulse of arbitrary 
form, the generated state \fu) — C/[(/)]F^|<I>f) can 
be represented as a coherent bosonic state as well, 
\fu) = nfe>o's"'^'~""^"^^l*B) with amplitudes Uk = 
vpVk 4'{kv-p)/2Tr determined by the Fourier transform 
(j){uj) of the phase (j){t) = f dt'eV{t')/h. The volt- 
age pulse V{t) generating the wave function behind 
the scatterer then derives from the phase 4>{uj) — 
(j){uj) exp[iSsc{uj/vF)]. 



S\fg) = exp(z$)J7i[x(i)]C/2[x'(i)]|/5) (42) 

with $ an overall phase. As a result, the scattered state 
is a Slater determinant and thus is non-entangled. 

In the following, we assume for simplicity that the in- 
teracting leads are characterized by the same coordinate 
kernels ki{x) — K2{x) — k{x). Let Ci, = {61,1,61/2} and 
e^, V = 1,2^ he the eigenvectors and corresponding eigen- 
values of the coupling matrix of the interaction Hamilto- 
nian Eq. (2). Provided the wave-packets f{x) and g{x) 
are Lorentzians with widths ^1 and ^2, the interaction 
equivalent phases x(t) arid x'{t) have Fourier components 
of the form, 
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X{uj > 0) = eii(eiie-'^^i + 6126"'^^=^ + e2i(e2ie-"^i + e22e"^^')F2{uj), 
X'{0J > 0) = ei2(eiie-""i + ei2e-""^)i^i(cj) + e22(e2ie-'""i + e22e-'""^)i^2(w), 

with x(— w) = X*('^) ^'^d the same for x'i'-^) ^i^d the overall phase reads 



did 

2^ 



lRe[G++{w)]\F,{u 



Re[F,{iu)] 



e,;ie 



+ 6^26" 



(43) 
(44) 



(45) 



where T1.2 = Ci,2/i'f- The functions Fi{u)) are given by 



F,{u) 



l-ie,n* (w)' 



(46) 



with k(cj) = / dt K{vpt)e^'^'^* , n++(u;) = j dxdx' k{x) 
k[x') G++{ijj,x,x') and G++(w, a;, a;') is the Fourier 
transform of the Green's function G_|_+(r, x, x') = 

{f+{p{x,T)p{x',m. 

G++iiu,x,x') = — ^ / — 7^-—-. (47) 



27rwp y 27rz w' — cj — i&gn(cj) 



C. The complete scattering state 

Due to the self-capacitance of the interferometer arms, 
electron-hole pairs can be excited in scattering processes 
where only one wave-packet (the upper or lower) propa- 
gates through the corresponding interacting arm of the 
interferometer. The corresponding interaction equivalent 
fields X and x' and the overall phase $ can be formally 
found from Eqs. (43)-(45) by setting the width of the 
absent wave-packet to infinity, ^1 — >■ cx) or ^2 ~^ 00. In 
the following, we introduce a new notation for the inter- 
action equivalent fields in the upper {xij ) and lower {xij ) 
interferometer and for the overall phase {^ij), where the 
indices i,j e {0,1} tell whether a particle has moved 
through the upper (i = 1) or lower (j = 1) lead (else 
i,j = if there is no particle in the lead). The scatter- 
ing state inside the interferometers before the outgoing 
beam splitters B and B' then takes the form (we assume 
symmetric incoming beam splitters with a = a' — 7r/4; 
note that the operators U always act on the interacting 
leads '1' and '2'), 



BB') - ^[M + e'''''U^[xn]U2[Xn]fl9l 



ft At 



Ui[xio]U2[x[o]fUl 



+e'*-f/i[X0i]C/2[x[,i]/l3^]l$F). (48) 

Each of the four scattered- or pseudo-spin components 
corresponds to a Slater determinant state and thus does 
not exhibit any interaction-induced entanglement. How- 
ever, the overall state cannot be factorized with respect 
to the upper and lower interferometer degrees of freedom 



because each component is formally obtained under a dif- 
ferent evolution (different interaction-equivalent phases). 
Hence, the interaction induces pseudo-spin entanglement 
between the non-entangled many-particle electron-hole 
cloud states. 

Before evaluating the Bell inequalities for the full setup 
including the Fermi sea, we first wish to gauge our expec- 
tations with a discussion of the degree of entanglement 
we can expect in our two Mach-Zehnder interferometers. 
Our Bell test is sensitive only to the excess particle num- 
ber transmitted through a given outgoing lead. In order 
to describe the statistics of the Bell test, it is sufficient 
to know the reduced 4x4 density matrix pbb' obtained 
from the full density operator = \BB'){BB'\ by trac- 
ing out all spatial degrees of freedom, see Eq. (6). For- 
mally, this corresponds to dividing all degrees of freedom 
into the (obervable) pseudo-spin degrees of freedom as- 
sociated with finding a particle in the leads '1' or '3' 
(encoded by B) or in the leads '2' or '4' (encoded in 
B') plus all the remaining degrees of freedom (we encode 
these with the letter R), hence pbb' = Ttr\BB') {BB'\. 
In a first step, we then can calculate the von Neumann 
entropy S{pbb') = -Tr[pBB/ log2(pB_B')] of the mea- 
sured outcomes. This entropy ranges between and 2 
and tells us about the degree of entanglement between 
the pseudo-spin degrees of freedom and the rest of the 
system. A value S[pbb') — tells us, that R is not en- 
tangled with the pseudo-spin system BB' and hence the 
system BB' can be maximally entangled; on the con- 
trary, a value S{pbb') ~ 2 informs us that R is fully 
entangled with BB' and all information about the latter 
is encoded in R. Due to the monogamy of entanglement, 
the two subsystems B and B' cannot be entangled in this 
case. Indeed, calculating the reduced densitiy matrices 



and the associ- 
(and similar for 



Pb = T^t:b'[pbb'] and ps' = T!t:b[pbb' 
ated entropy S{pb) = -Tr[/3B \og2{pB) 
B') we expect that S{pb) ranges between {B not en- 
tangled with B') and 1 (B maximally entangled with B') 
if S{pBB') = 0, while S{pb) = if S{pbb') = 2. 

Hence, the more degrees of freedom from R that are 
entangled with B or B' we integrate out (as we do not 
observe them in our Bell test), the smaller is the remain- 
ing entanglement between B and B' and the weaker is 
the violation of the Bell test, see Fig. 4. Consider first 
the two-particle Bell test described in Sec. III. For nar- 
row wave packets, the trace over the coordinates does 
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FIG. 4: Von Neumann entropy S{pbb') of measured out- 
comes of the double Mach-Zehnder interferometer quantify- 
ing the entanglement with the unobserved degrees of freedom 
in the system. The entanglement increases with the inter- 
action strength and the system is fully entangled with the 
environment at large ifio- The entropy S{pb) expresses the 
entanglement between the two Mach-Zehnder interferometers 
and thus the potential for violation of the Bell test. Increas- 
ing the interaction first entangles the two MZ interferometers 
with each other but as the entire system gets entangled with 
the environment upon further increase of (po, the entangle- 
ment between the interferometers decreases. 



not reduce the entanglement in the outgoing leads (hence 
S{pbb') — 0) and we expect the entropy S{pb) to reach 
unity under ideal conditions where B and B' are fully 
entangled; we then can expect that the Bell inequality 
can be maximally violated. For extended wave packets, 
the spacial degrees of freedom are relevant, e.g., some 
'which-path' information is stored in the shape of the 
wave-function in the outgoing channels. The entropy 
S(pbb') does not vanish as the system BB' is entan- 
gled with R and hence the entropy S{pb) cannot reach 
unity any more; the Bell inequality cannot be maximally 
violated any longer. Finally, including the Fermi sea, the 
entanglement between R and BB' is even larger since 
even further information on the BB' system is encoded 
in unobserved degrees of freedom in R, e.g., the pres- 
ence of an additional hole in the outgoing channel tells, 
that the particles sure did not choose the paths '3' and 
'4', since no interaction is active in this case. As a conse- 
quence, S{pbb') deviates more strongly from zero, S(pb) 
is further reduced from unity, and the maximum in the 
Bell inequality violation is further diminished. 

For a quantitative analysis, we then calculate the 
two-particle density matrix pbb' corresponding to the 
state (48), see Eq. (6). The diagonal elements are all 
equal to 1/4, while the remaining six independent ele- 
ments are given by 



[/5bB']i2,14 
[PSB']l2.32 



-1 f T2dt ^^[x;l]+(^)+^[x'lol-(*) 



Axil]+(t)+i[X0l]-(t) 



47r J ri+t^ 
-I f Tidt 



47r J rf+t^ 



(49) 



(50) 



[PBB'] 14,34 



--cxp(i[xio] + («Ti)), 

-^CXp(i[xoi] + («T2)), 



jnrf^ ,[xio] + (t)+*[xoi]-(t) 



T2dt 



(52) 
(53) 
(54) 



rl 



-«[Xoi]-(*)-«[xio] + (*) 



where [/]±(t) are the components of the function f{t) 
analytical in the upper (lower) half plane. 



1 



dt' 



fit') 



t' -tTiO 



(55) 



For the sake of simplicity we assume a symmetric sit- 
uation where the interacting leads have the same self- 
interaction Ell — E22 and the wave-packets have the 
same form f{x) ~ g{x) and thus ti = T2 = r, im- 
plying that [pBB']i2,ii = [PBS']i2,32 ^nd [pBB']i4,34 = 

[ps_B']32 34- We first focus on the situation of infinitely 
narrow wave-packets 7 = ^/a — )■ 0, allowing us to re- 
place r/(r2 + e) TrS{t) in Eqs. (49)-(54). The Bell 
correlation function then takes the form 



E = 



l2/3cos2/3'[cos($ + (/7s+<Pc) - cos(<i> + (^s)] 



^ cos2/?sin2/3'[cos($' + (p5 + (^c) - cos($' + </?s)] 

- isin2/3sin2/3'[(ViV^2)^cos($ + $' + 2(/:>, + (^e) 

+ (V^i/X^2)'cos($-$')]- (56) 

Here, the phase (pc = Re[xoi] + (^ = 0) assumes the role of 
ipc in the two-particle case, the phase difference between 
the two cases where a particle traverses the lead '1' with 
(Xii) and without (xio) presence of a particle in lead 
'2'. Since Xii — Xio + Xoij it is the difference xoi that 
enters into ipc- The phase ips = Re[xio]-i-(i = 0) is due to 
the self-interaction of the particle passing through lead 
'1'; this phase is absent in the two-particle scattering 
problem. Finally, the visibility factors are given by Vi — 
exp[— Im[xio] + (i = 0)] (describing the situation where 
only one particle traverses the interaction region) and 
V2 = exp[— Im[xoi] + (i = 0)] (showing up when both 
particles traverse the interaction region; no correlations 
appear in the absence of any interaction). Note that both 
factors ^1^2 < 1 and V1/V2 < 1. 

The phase shift (ps can be eliminated from Eq. (56) by 
incorporation in the magnetic fluxes and setting ips — 0. 
Let us then choose the same values of (renormalized) 
magnetic fluxes as for the two-particle case, Eq. (29); 
substituting the settings in Eq. (29) into Eq. (56), we 
find the Bell parameter 




sm 



Pc 



(57) 



[PBB' 



at ^ = 7r/4 and cot 2/3' = 2 sin((^c/2)/(Fi(y2^ + V^2~^))- 



12,34 



- expii[Xii] + {iTi) + i[x'ii] + iiT2)) , (51) ^s before, we consider the interaction kernel k{x) 
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FIG. 5: Bell parameter (57) for the non-optimal settings of 
the Bell parameters in Eq. (29) as a function of interaction 
parameter (po (horizontal axis) and the parameter rj (vertical 
axis) describing the ratio of self- to mutual capacitance. 



FIG. 6: Maximal value of the Bell parameter as a function of 
interaction parameter (po (horizontal axis) and the parameter 
•q describing the ratio of self- to mutual capacitance (vertical 
axis) . 



exp(— |a;|/a) and relate the self-coupling constants En to 
the mutual coupling parameter E12 via the dimensionless 
factor ?7 > 0, £'11 = E22 



and 



rjEi2. Then 
1 + {uJTaY ' 



To 



27ri (wtq + lY 



(58) 



(59) 



with To — a/vp the ballistic traveling time through the 
interaction region. For the symmetric setup, the inter- 
action equivalent phases obey the symmetry relations 
Xio = Xoi and xoi = x'lo '^ith 

/^-\u\i+iiyt/To 
di^-, — —^^77 (60) 

where nio = (i^o/27r)(7]+l), and similarly for the xoi field 
with noi = {(po/2TT){ri — 1). In Fig. 5 we plot the value 
of the Bell parameter (57) as a function of ipo, see Eq. 
(22), and the parameter 77. The maximal value B « 2.13 
violating the Bell inequality is assumed at tpo ~ 0.53 tt 
and r] = 0.083. 

As noticed before, the settings (29) might be non- 
optimal in the general case. The full optimization of the 
Bell test according to Eq. (7) gives a higher value of the 
Bell parameter, see Fig. 6, where we find the maximal vi- 
olation of the Bell inequality B w 2.18 > 2 at (ySo ~ 0.73 tt 
and 77 « 0.58. Note that a value 77 < 1 implies that the 
self-capacitance is smaller than the mutual capacitance, 
E22 < E12; we will discuss this point further below. 

Above, we have concentrated on the case of infinitely 
narrow (or high energy) wave-packets, 7 = ^/a — > 0. 
The maximal value of the Bell parameter for the general 
situation with a finite- width wave-packet is shown in Fig. 



7, where we have chosen the optimal value rj = 0.58 for 
narrow wave packets with 7 = 0. We find that lowering 
the energy of the excess electrons decreases the value of 
the Bell parameter and using high-energy wave-packets 
happens to be the best operating limit for the Bell test. 

The interaction-equivalent phase xw involves two di- 
mensionless parameters, 7 = ^/0 and the Coulomb 
strength ipo = Ei2To/h. Using phase-space arguments, 
one can see that these parameters determine the num- 
ber of excitations that can be created via the wave- 
packet's excess energy de^ = Hvp/^S^ or via the Coulomb 
energy Sec = £'12 • Indeed, multiplying these energies 
by the density of states p = 1/hvp and the length a, 
we obtain the number of electrons that can be excited 
within the interaction region, ^ apSe^ — a/2^ and 
Nq — SecT^/h = ipo- The narrower the wave-packet 
is, the more electron-hole pairs can be excited. On the 
other hand, the finite Coulomb interaction strength re- 
stricts the number of excited particles to below iVc, no 
matter how narrow the incoming wave-packet is. As a 
result a further increase in the energy of the incoming 
particles beyond i5ec is not harmful since it does not pro- 
duce more electron-hole excitations and thus does not 
lead to a further reduction of the Bell parameter. 

Decreasing the number of parasitic electron-hole pairs 
should be a good strategy to enhance the value of the 
Bell parameter. Keeping the interaction fixed (say, near 
ipa = tt), one could choose spatially extended wave- 
packets such that iV^ drops below Nq- In this situation 
the number of electron-hole pairs is defined by the smaller 
number N^. However, for extended wave-packets with 
large f , the electrons are unlikely to meet each other in 
the interaction region, thereby decreasing the interaction 
phase, see Eq. (21). As a result of this trade-off, the Bell 
inequality can never be violated in this regime, though 
the number of parasitic electron-hole excitations is small. 
However, using extended wave-packets, one could count 
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FIG. 7: Bell parameter (57) as a function of interaction pa- 
rameter ipo (horizontal axis) and the width parameter 7 = ^/0 
(vertical axis) at 77 = 0.58. 



FIG. 8: Interaction region between edge states '1' and '2' in 
a quantum Hall device. The edge states are separated from 
the closest screening Fermi sea by a distance i'^ (top 
left). We model the capacitive interaction Hamiltonian (2) by 
analyzing two charges qi and 52 separated by s and screened 
by metallic plates at a distance dp. 



the particles in the outgoing leads within a finite time 
window St, Ni{t) — J^^g^j^ dt ii{t), thus projecting the 
electron trajectories on the component where the elec- 
trons have passed the interaction region simultaneously. 
For ^ > a, the number of excited electron- hole pairs is 
limited by <C 1 and one can observe a sufficient vio- 
lation of the Bell inequality. In this approach, however, 
the entanglement is likely to be induced by the projection 
(or the Bell measurement) itself and the setup cannot be 
used as a genuine source of entangled particles. 

So far, our results have been obtained for the specific 
form of the interaction kernel k{x) = exp(— |x|/a) with 
smooth tails. A different form of k{x) may result in a 
different value of the maximal Bell parameter. To check 
how strongly the result depends on the shape of k(x), we 
consider another kernel k{x) = 1 for x G [—a, a] and 
otherwise. This kernel has sharp edges and generates an 
interaction equivalent phase of the form 



component 



Xw{t) = 4nio / diy 



2 

sm V 



-\iy\l+ii^t/To 



V V ~ mio(l — e '^™) 



(61) 



and similarly for the phase xoi with tiiq — > ^oi, see Eq. 
(60). The numerical analysis for high-energy wave pack- 
ets (7 — )■ 0) produces a maximal value of the Bell pa- 
rameter B ~ 2.15 which is assumed at lower interaction 
<y9o ~ 0.34 TT and a similar value r\ ~ 0.53. 

Finally, we discuss the capacitive Hamiltonian (2) used 
throughout the discussion. In a realistic situation the in- 
jected electrons experience a Coulomb interaction which 
is screened by the Fermi sea when propagating through 
the adjacent interaction leads. The peculiarity of the 
edge states of the Quantum Hall system is that the edge 
state ipk{y) e''^^ with a wave- vector A: > has a transverse 



(.) = ^«(^)exp( 



{y + 



(62) 



which is spatially separated by from the Fermi sea 
electrons occupying states with fc < 0, see Fig. 8 (here 
£ — \J ip^l^TiB is the cyclotron length and Hn{y) are Her- 
mite polynomials; we assume a flat boundary potential). 
Denoting by dp the distance between the two Fermi seas, 
Lorentzian wave-packets with a small width ^ then can 
be brought close to each other, with a distance dp — 2ds 
between the two wave packets and a distance d^ « ^^/C 
from the nearest Fermi surface screening the interaction. 

We model this situation by two parallel grounded 
metallic gates (accounting for the screening effect of the 
Fermi seas in the adjacent interacting arms) separated by 
a distance dp with two charges qi and 52 at a separation 
s < dp and located symmetrically along a line perpen- 
dicular to the metallic plates. The electrostatic energy 
of this system is given by, 



E{qi,q2) = + ^12 9i'Z2 

with coupling constants (£^11 = £^22) 



£-22 2 
^92, 



(63) 



Eii 
E12 



1 

dp 
1 

dp 



,(i±£),,(i^) 



(64) 



where ^'{x) is the digamma function and a = s/dp. The 
last term cx 1/a in the off-diagonal coupling strength 
E12 is the direct Coulomb interaction qiq2/s of the two 
charges. It is this term which allows to realize a situation 
where |£^i2| > \Eii\ when a < 0.177. 
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V. CONCLUSION 

We have studied the on-demand generation of entan- 
gled electron pairs in a double Mach-Zehnder interferom- 
eter implemented in a quantum Hall setting; the entan- 
glement is generated via the capacitive interaction be- 
tween two neighboring interferometer arms, shifting the 
phase of one wave-function component with respect to 
the other. The resulting entangled state is analyzed in a 
Bell test measuring the particle numbers in the outgoing 
leads; the second beam splitters and the fluxes thread- 
ing the interferometers serve to define the four measuring 
conditions (polarizations) in the Bell test. 

Our first study with only two electrons present in the 
device serves to identify the optimal conditions to gener- 
ate full entanglement and maximal violation of the Bell 
test in the simplest situation. We find that simultaneous 
injection of high-energy pairs provides the best condi- 
tions for entanglement; broad (low-energy) wave-packets 
reduce the probability that the electrons meet and inter- 
act, as does a delay between the particles. Furthermore, 
extended wave-functions are deformed in the interaction 
region and thus transfer 'which-path' information to the 
outgoing leads. Since our Bell test is sensitive only to the 
particle number in the outgoing leads, this (discarded) 
wave-function deformation entails a reduction of viola- 
tion in the Bell test. 

Including the Fermi sea, decoherence and addi- 
tional 'which-path' information due to the generation of 
electron-hole pairs within the interaction region are ex- 
pected to further reduce the entanglement in the relevant 
degrees of freedom, the particle numbers measured in the 
outgoing leads. The most important question addressed 
in the present work then is: can such an idealized de- 
vice generate sufficient entanglement to be observed in a 
Bell test ? We find that this is indeed the case, provided 
the capacitive coupling between the two arms can be im- 
plemented such that the mutual interaction dominantes 
over the self-interaction. 

In our work we consider electrons in chiral states of 
a quantum Hall device and a purely capacitive (in par- 
ticular, non-resistive) coupling. Technically, this allows 
for the bosonization of the scattering problem describ- 
ing the interaction region (however, not the scattering at 



the beam splitters) , which is at the heart of transforming 
an incoming Slater determinant into an outgoing Slater 
determinant; this allows us to describe the effect of the 
interaction via a voltage-pulse, thus replacing an appar- 
ent many-body problem by a single-particle evolution. 
Physically, this implies that the decoherence is reduced 
in our system and limited to the creation of particle-hole 
pairs. It turns out, that these particle-hole pairs further 
reduce (as compared with the two-particles case) the en- 
tanglement and the violation in the Bell test, but to a 
limited degree, still providing a violation B w 2.18 > 2 in 
an optimized situation. Regarding the experimental im- 
plementation, we refer the reader to the layout described 
in Ref. Iti. A favorable finding is that the Bell test can 
be performed by only changing the fluxes through the 
loops (possibly by using a side gate deforming the loop 
area) and does not require tuning of the beam splitters. 
An interesting further element is, how correction pulses 
resurrecting the wave function behind the scatterer as 
described in Ref. 22 can help to restore a higher degree 
of Bell-inequality violation. 

We acknowledge financial support from the Swiss Na- 
tional Science Foundation through the National Center of 
Competence in Research on Quantum Science and Tech- 
nology (QSIT), the Pauli Center for Theoretical Studies 
at ETH Zurich, and the RFBR Grant No. 11-02-00744-a. 



Appendix A: Wave packet scattering 

We show that the scattering state of two Lorentzian 
wave-packets fi^2{x) propagating along two interacting 
adjacent leads '1' and '2' is indeed described by a Slater 
determinant state and thus can be prepared by apply- 
ing an evolution operator generated by a single-particle 
Hamiltonian. We assume that this single-particle evolu- 
tion can be generated by voltage pulses, i.e., the Hamil- 
tonian accounts for the interaction of the electrons with 
voltage pulses applied at x = and adds the time depen- 
dent phases x(t) and x'(t) to the state in the leads '1' and 
'2', see Eq. (34). In order to demonstrate that there are 
specific phases x(^) a-^d x'(0 generating the scattering 
state, we need to calculate the overlap. 



= j dxidx2rAxi)f2ix2) J dx[dxy,{x[)f2{x'^){^2{x2,to)i'iixi,to)ul[x]uUx']S^^^^ (Al) 

and show, that there exist phases xi^) and x'{t) such that the modulus \0\ reaches unity. In Eq. (Al), to — )■ — oo and 
the averaging is taking over zero temperature Fermi sea. In order to calculate the overlap, we can make use of the 
Green's function, 

C(12|1'2') - {TK{Ui[x{t-)]U2[x{t-)]S+^2{x2,to)^i{xi,t^M{x[,4M^^^ (A2) 

I 

where we have introduced the Keldysh time ordering Tk branch of the Keldysh contour (going back in time) , while 
and have defined the fields x(t) and x'(i) on the lower 
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the evolution operator 

S+ = Tk exp 



+ 00 



dtHintit+) 



(A3) 



is defined on the upper branch of the Keldysh contour 
(forward in time). 

We decouple the quadratic interaction Hamiltonian (2) 
with the Hubbard-Stratonovich transform Eq. (13) and 
obtain the evolution operator S+ in the form 



^4 



(A4) 



X exp 



dt[eizl{t) ^- t2zl{t)] 



where the real fields zi.2(t) are non-zero on the upper 
branch of the Keldysh contour. The actions (see Eq. 
(12) for the definition of rii) 



Si[zi\ = Tftrexp -i / dt Zi{t)n.i{t) 



1,2, 



(A5) 



with Zi(t) = eiZi{t)en + e2^2(i)e2i, describe the forward 
in time evolution of the electrons in the leads i = 1, 2 
under the action of the time dependent fields Zi{t) acting 



within the interaction region according to the Hamiltoni- 
ans Hzi(t) = J dx Zi{t)K{x)pi{x). The Green's function 
(A2) can be factorized with respect to the leads, 

C(12|1'2')-(C(1|1')C(2|2')),^_,^, (A6) 

where the average is taken over the fields zi 2 and 

CixtoWto) = {Tk{U[x-]S[?+] (A7) 
xi'{x,t^)¥{x',t+)}) 

with S'[z+] the action in the relevant lead. Next, we 
make use of the bosonization technique, see Ref. 2') for 
details, and express the density operator p(x, t) through 
the chiral bosonic field p{x, t) = dx9{x — vpt) /27r where 

Hx) = -Efc>o(^fce'''" +ble-'''^)/Vk and bl bk are 
bosonic creation and annihilation operators obeying the 
standard commutation relation [bk, b\,] ~ 5kk' ■ The elec- 
tronic field operator '^{x) can be expressed through the 
field 0{x) via 4'(a;) = Fe'^^^^Vv^ with 5 ^ 0+ an 
ultraviolet cutoff and F is a Klein factor acting as an 
electron-number ladder operator. 

Calculating the bosonic averages in (A7), one arrives 
at the result 



C(r|r') 



1 (>i[x\+{r')+i[x]~{r) 

2TTivF T — t' — i6 



■ exp 



duj r 

2^ 



|z(c.)|2n++(^)-z;i|x(c.)pG;+(c.) 



where we have introduced retarded variables t = to — x/v-p and t' = Iq — x' /v-p; the quantities G^+{uj) and Il++{uj) 
are defined in Eq. (47). Substituting this expression into Eq. (A6) and taking the Gaussian integrals over the fields 
zi{t) and Z2(t), one finally arrives at the following expression for the correlation function C(12|1'2'), 



e'-lx]+{T[)+i[x]-(ri)+i[x']+{ri)+i[x']-ir2) 

C(12|1'2') = _, ... exp 



X exp 

fc=1.2 



(27riwF)^('''i ~ 't'i^ *'5)(t2 — ^2" iS) 
duj 
2 



duj 
2^ 



(|x(a.)P + |x'(c^)P)G;+(a.) 



(A9) 



with the functions Fj.(cj), k = 1,2, defined in Eq. (46). 
This correlator is an analytic function of the retarded 
variables ri and T2 (t( and Tj) in the lower (upper) half 
plane. This feature allows us to evaluate the integrals 



over the coordinates in Eq. (Al) involving Lorentzian 
wave-packets with the widths ^1 and ^2- In the end, we 
arrive at the following expression for the overlap 
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O = exp 



I rpi\x{-)\'+\x'{-^)\')GU{^) 



X exp 

fe=l,2 



2^ 



Minimizing the real part of the overall exponential pro- producing a maximal overlap \0\ = 1; the overall phase 
vides the optimal fields x and x' in Eqs. (43) and (44) is given by Eq. (45). 
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